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I. INTRODUCTION 

A growing interest has recently been devoted to 
the study of quantum phase transitions in spin 
systems especially from the entanglement point of 
view. The ground-state intricatcd structure has in- 
deed been shown to be strongly affected by the ex- 
istence of a critical point as initially illustrated in 
the one-dimensional (ID) Ising model under magnetic 
fieldi^. Following these pioneering works, many sys- 
tems have been studied and have revealed a specific 
behavior of the entanglement as measured either by 
the Concurrenco3,4,5,6,7,8,9,io,ii,i2,i3,i4,i5,i6,i7 by the 

entropy^^i^^i^°i^^^'^i'^'^i^'^. In this context, the Lipkin- 
Meshkov-Glick (LMG) modelSSiSSiSi has focused much 
attention because of its apparent simplicity. Introduced 
fourty years ago in nuclear physics, this model provides 
a simple description of the tunneling of bosons between 
two degenerate levels and can thus be used to describe 
many physical systems such as two-mode Bose-Einstein 
condensates^^ or Joscphson junctions. In the thermody- 
namical limit, its phase diagram can be simply estab- 
lished using a semiclassical (mean-field) approach22iffi. 
However, at finite size, the situation is more complicated 
and crucial to understand the entanglement properties of 
this niodeli^iSi^SiSi^. A possible way to investigate this 
problem is to use the exact solution provided by the al- 
gebraic Bethe ansata^^^, but it would require the com- 
putation of the correlation functions which is a rather 
difficult task with this formalism. 

In the present work, we present an alternative route 
which relies on a combination of two well-known meth- 
ods, and which we have already briefly sketched in 
Ref. 13 First, we use the Holstein-Primakoff represen- 
tation which allows us to write a 1/N expansion of the 
spin operators at arbitrary order, N being the number of 
degrees of freedom. Usually, such a development is often 
restricted to the first order for which the Hamiltonian re- 
mains quadratic and can thus be easily diagonalized. We 
show that this first order calculation does not provide us 



with enough information. This implies one has to take 
more terms in the 1 /N expansion, and thus to deal with a 
nonquadratic Hamiltonian. Here, we use the continuous 
unitary transformations (CUTs) technique proposed by 
Wegner— and independently by Glazek and Wilson^S*^, 
to compute the next orders. The structure of this ex- 
pansion contains all relevant informations to extract the 
finite size corrections of various obscrvablcs at the crit- 
ical point. We consider in particular the ground-state 
energy, the gap, the magnetization, and the spin-spin 
correlation function for which we obtain nontrivial scal- 
ing exponents. These results corroborate several numer- 
ical studiesi^iSS in which some of these exponents have 
already been computed. 

In addition to giving details of the calculations yielding 
the results announced in Ref. l34l we present a numeri- 
cal study for all the finite-size scaling exponents which 
corroborates the analytical predictions. Furthermore, we 
investigate the broken phase whereas in Ref. |3 we only 
considered the symmetric one. Finally, we compute the 
leading order two-spin entanglement in the thermody- 
namical limit for any anisotropy parameter. 

This paper is organized as follows. In Sec. ^ we in- 
troduce the LMG model, discuss its symmetries and give 
its simple solution for the isotropic case. In Sec. IIIII we 
recall the variational approach which constitutes the ze- 
roth order of the 1/N expansion but which is sufficient 
to establish the phase diagram in the thermodynamical 
limit. Section Hvl is devoted to the first order corrections 
that are easily computed via the standard Bogoliubov 
transformation. In the next section (Sec. we intro- 
duce the basics of the CUT formalism and apply it to 
the LMG Hamiltonian. These corrections allow us to 
give analytical expressions of the finite-size scaling ex- 
ponents. In Sec. IVII we present a numerical study that 
confirms our results and sheds light on the discrepancy 
with previous numerics. Then, in Sec. lVIll we discuss the 
entanglement properties of the ground state focusing on 
the so-called concurrence which measures the two-spin 
entanglement. We show, in particular, that it displays a 
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cusp at the transition point contrary to what is observed 
in 11? systems. Technical details are given in Appendix. 



II. THE LIPKIN-MESHKOV-GLICK MODEL 



A. The Hamiltonian and its symmetries 

The LMG modctfSiSSiSI describes a set of spins half 
mutually interacting in the (anisotropic) x ~ y plane em- 
bedded in a perpendicular magnetic field pointing in the 
z direction. The corresponding Hamiltonian is written 



i<3 



2A 
A 



~-^{Sl + lSl)-2hS, + -{l+^l (2) 



--(1 + 7) {S^-Sl-N/2)-2hS, 



(3) 



where the ctq's arc the Pauli matrices, Sa = Tlii'^al'^^ 
and S± = Sx ± iSy . The prefactor 1 /N is necessary to 
obtain a finite free energy per spin in the thermodynam- 
ical limit. For any anisotropy parameter 7, H preserves 
the magnitude of the total spin and does not couple states 
having a different number of spins pointing in the field 
direction (spin-flip symmetry), namely 



[H, S^] = and 



(4) 



An important consequence of the spin-flip symmetry 
is that any cigcnstatcs of Yli satisfy 



{Sx) = {Sy) = 



and 



{SxSz 

{SySz 



(SzSx) 

(SzSy) 



(5) 



(6) 
(7) 



In the following, we only consider a ferromagnetic dom- 
inant coupling (A > 0, 17I < 1) and without loss of gener- 
ality, we set A = 1. For a discussion of the antiferromag- 
netic case, we refer the reader to Ref.ls^ We also restrict 
our discussion to the region h > since the spectrum is 
even under the transformation h —h. In addition, we 
only consider the maximum spin sector S ~ N/2 which 
contains the low-energy states. 

The ground-state properties can be easily studied in 
the thermodynamical limit using a mean-field analysis 
that we briefiy recall in the following section. It is worth 
however to first have a look at the isotropic model for 
which a complete analytical solution is available. 
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FIG. 1: The mean-field ground state in the symmetric phase, 
for all values of the anisotropy parameter 7. 



B. Analytical solution of the isotropic model 

In the isotropic case 7=1, the Hamiltonian ^ reads 

A^' 



H 



2 



2hSz. 



(8) 



H thus commutes both with and so that it is 
diagonal in the standard eigenbasis {\S,M)} of and 
Sz- The cigenenergies are 



EiS,M) = --[S{S+l) 



- N/2] - 2hM. (9) 



For the ferromagnetic case that we consider here, the 
ground state is then obtained for S = N/2 and 



I{hN/2) for < /i< 1 
N/2 for h > 1 



(10) 



where I{x) gives the integer part of a;, in the following 
sense : ii x ^ X + 5x with X an integer and (5a; G [0, 1[, 
then I{x) = X for Sx G [0, l/2[ and I{x) = X + 1 for 
6x G [1/2, 1[. One then trivially gets 
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(16) 



where eo(N) is the ground-state energy per spin. The 
gap can also be easily computed : 

A(Af) = E{N/2, Mo + -s) - E{N/2, A/q) (17) 

= ^ + 2s{2Mo/N -h), (18) 
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FIG. 2: The two mean-field ground states in the broken phase 
of the anisotropic LMG model 7 < 1. 



FIG. 3: One of the infinitely degenerate mean-field ground 
states in the broken phase of the isotropic LMG model 7 = 1. 



where s = sgn{hN/2 ~ Mq) ii < h < 1 and s = -1 if 
h>l. This gives A{N) = 2/N + 2s{2Mo/N - h) in the 
broken phase [which is a 2 /A^-periodic function vanishing 
for h = {2p+ l)/N and equal to 2/N for h = 2p/N, p 
denoting an integer], whereas A (TV) = 2{h — 1) + 2 / N in 
the symmetric phase. 



III. VARIATIONAL AND SEMICLASSICAL 
APPROACHES 

In this section, we describe the semiclassical approach 
that can be used to determined the phase diagram of 
the LMG modeli22i2£ This analysis constitutes the zeroth 
order approximation and relies on a mean-field (varia- 
tional) wave function 

N 

\m </>)) - (g) [cos (9/2) e-'^/2| t>,+sin (0/2) e'*/'\i)i] , 
1=1 

(19) 

which is a coherent spin state such that 
N 

(S) = — (sin0 cos 0, sin6' sin (/), cos 0). (20) 

In Eq. ((T^ . kets 1 1)/ and | j); are the eigenstates of (t[ 
with eigenvalues -fl and —1, respectively. The ground 
state is thus determined by minimizing the energy 

(H) = - sin2 e (cos^ (/) + 7 sin^ 0) - hN cos 0. 

^ (21) 
with respect to 9 and 0. 

In the thermodynamical limit, this variational ap- 
proach is completely equivalent to the semiclassical treat- 
ment proposed by Botet and JuUien^ which consists in 
treating S as a classical variable. In this limit, minimiz- 
ing {H) with respect to 9 and (j) leads to a distinction 
between the two phases 

1. h > 1 (symmetric phase): The ground state is 
unique and fully polarized in the magnetic field di- 
rection {9o = 0) for all 7, see Fig.^ 



2. Q < h < 1 (broken phase): 

For 7 7^ 1, the ground state is twofold degenerate 
{9o = arccos/i and 4>q = or n) (see Fig. EJ. In 
the isotropic case 7 = 1, (H) does not depend on 
(j) so that the ground state is infinitely degenerate, 
as shown in Fig. |31 Of course, this is due to the 
parametrization of the trial function. However, if 
this result is clearly in contradiction with the exact 
degeneracy computed in Sec. IIIBI (which is unique 
at least for /i 7^ 0), it points out the existence of 
two universality classes (7=1 and 7 7^ 1). 

Here, 9o and (f>Q stands for the values of 9 and for which 
(H) is minimum. The ground-state energy per spin is 
then given by: 



Co (00) = — - sin^ 9o — /i cos 6*0, 



(22) 



and the gap A can be evaluated using the random phase 
approximation^ 



A(oo) 



^ 2[ih-l)ih-j] 
I 



,1/2 



for h>l, 
for 0<h<l. 



(23) 



Let us note that, for the anisotropic LMG model, the 
result in the broken phase < h < 1 simply says that 
there are two degenerate ground states, and is obtained 
without performing a calculation. This gap does not cor- 
respond to the excitation energy around one of the two 
mean-field states represented in Fig. |21 which is finite. 
The latter is denoted as A' in Fig. ^ whereas in the 
same figure, A denotes the exact gap, vanishing in the 
broken phase and in the thermodynamical limit. 

The correlation functions are also easily obtained in 
the thermodynamical limit since, for a "classical" spin 
one simply has {ScSff) = (5'«)(5'/3) + 0{N). Note that 
for 7 = 1, one can recover the exact expressions derived 
in 111 HI provided an average over all possible values of 4>o 
is performed. For example, one has: 
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FIG. 4: Schematical representation of the low-energy spec- 
trum of the anisotropic LMG model in a finite system. A is 
the energy gap between the ground state and the first excited 
state. In the broken phase h < 1, A is exponentially small 
A ^ exp(—aN)^ In this phase, A' is the excitation energy 
between the two nearly degerate lowest states and the nearly 
degenerate next two states. 



n'Oo cos2 0o + 0(l/iV), 
■sin^ 6*0 + 0(l/iV). 



(25) 
(26) 



which coincides with Eq. H15() in the infinite N limit. 

In the thermodynamical limit, this mean-field approx- 
imation is exact (see the comparison with the numerical 
results in Sec. lVip and thus predicts a second-order quan- 
tum phase transition at the critical field h ~ 1. 

For finite N, the exact ground state is nondegenerate 
and is not given by a mean-field state {ipiS, (p)) except for 
h — -^^t^^/t- twofold degeneracy (for 7 7^ 1) in the 
broken phase only occurs in the infinite N limit where the 
two mean-field solutions \ip{Oo, 0)) and \tl.' {60,11)) provides 
a (nonorthogonal) basis of the ground-state subspace. 

Finally, from the entanglement point of view, 
\'ip{9o,(j)o)) is a completely separable state but as dis- 
cussed in Sec. lVIII the exact ground state has some non- 
trivial intrication properties (even in the thermodynam- 
ical limit) that cannot be captured at this level. It is 
thus essential to go beyond this basic approach and to 
compute next orders in the expansion. 



IV. FIRST ORDER QUANTUM CORRECTIONS 

We now wish to go one step beyond the mean-field 
analysis of the previous section by computing first order 
correction in a 1/iV expansion. To achieve this, we use 
the Holstein-Primakoff representation of the spin oper- 
ator and truncate the resulting bosonic Hamiltonian to 
lowest order. Next, we diagonalize it thanks to a Bogoli- 
ubov transformation. 



A. The Holstein-Primakoff representation 



The Holstein-Primakoff representation of spin 
operators^ is a useful tool to compute low-energy 



corrections around a classical magnetization (see for 
example Ref. |4^ . The very first thing to do is thus 
to perform a rotation of the spin operators around the 
y axis, that brings the z axis along the semiclassical 
magnetization. This is done as follows 



Sx\ ( cos 00 sin 6*0 

Sy\=\ 10 

Sz I V — sin 6*0 cos 6*0 




(27) 



As explained in the previous section, = for h > 1 
so that S = S, and = arccos/i for h <1. Note that 
we have chosen here the magnetization direction corre- 
sponding to 00 = but the same results can be obtained 
by choosing 0o = tt. 

The Hamiltonian written in terms of the S then 
reads 



H 



1+7 
2 
2 



2hmSz 



(1 - m'ysi + 



- Si 



-h\/l-m^ (S+ + S 



(28) 



N 



^S-\-Sz SzS-\- S—Sz ~\~ SzS— 



where m = 2{Sz)/N = cos^o- The Holstein-Primakoff 
representation is then applied to the rotated spin opera- 
tors 



Sz 



S-a''a = N/2 - a'' a, 



(29) 



S+ = (2S'-ata)i/2^^^i/2^^_^ta/7v)'^^a, (30) 



1/2 



S- 



t {2S - a^af^ = iV^/^at (1 _ a^a/NY'" , (31) 



1/2 



where the bosonic creation and annihilation operators 
satisfy [a, a^] = 1. 

The next step consists in inserting these expressions 
in Eq. H28|l . and to expand the argument of the square 
roots. Keeping terms of order (l/iV)~^, (l/iV)"^/^ and 
(l/iV)*^ in the Hamiltonian [therefore, it is sufficient to 
use the approximation (1 — a^a/NY^^ ~ 1] yields 

=iVeo(oo) + <5e[,^) + A(B)ata + r(B)(at' + a'), (32) 

where eo(oo) is the mean-field ground-state energy per 
spin (|22|) . which also reads eo(oo) = { — l — 2hm + m'^)/2. 
The superscript stands for "bare" , and the bare cou- 
plings are 



6e 



(B) 



9 ' 

A(^) = 2 + 2/iTO - 3™^ _ 



r(B) 



7 — m 



(33) 
(34) 

(35) 
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We emphasize that the Hamiltonian contains no term 
proportional to \^{a'^ + a), i.e., of order (l/iV)~^/^, cre- 
ating one excitation. These terms simply cancels because 
the angle of the rotation H27|) has been chosen to bring 
Sz along the classical magnetization. Note that this can 
also be achieved, as in Ref. 0, for example, by applying 
directly the Holstein-Primakoff mapping to the nonro- 
tated S operator, and performing a shift of and a by 
an appropriate constant. 



B. The Bogoliubov transformation 

The quadratic Hamiltonian H32|) is easily diagonalized 
by a standard Bogoliubov transformation 

at = cosh(e/2)fet + sinh(e/2)6, (36) 
a = sinh(e/2)5t + cosh(e/2)6. (37) 

where G is such that the Hamiltonian expressed in terms 
of the 6's does not contain a term 6^ -\-f) ^ i.e. is diagonal. 
Let us set e = -2r(^V^^^^ which takes the following 
values in the symmetric and broken phases 



e{h > 1) = 
e{0<h<l) = 



1-7 



2/1-1-7' 



■ 7 



2-/72-7 



(38) 
(39) 



One easily finds that to diagonalize the Hamiltonian H32|) 
one must satisfy tanh Q — s. One then has 



= iVeo(oo) + Se^^^ + A^^^b^b, 



(40) 



where the superscript means "renormalized" and the 
renormalized couplings are given by 



(B)^ (TT^-l), (41) 



A(R) 



.(B) 



(42) 



At this stage, it is worth noting that, in the broken 
phase, for h = y/^, e and thus F^^) vanishes so that 
the Hamiltonian H32(l is readily diagonal. This special 



point coincides, at this order (1/A^)°, with h 



jV-l 
N 



V7 



previously mentioned (see Sec. lHI|) for which the product 
state {ipidoj 4>o)) is the exact ground state. 
In the symmetric phase m = 1, one thus gets 



6e, 



(R) 



-/.+ i±^ + [(/.-l)(/^-7)]'/^ (43) 



A(^) = 2[(/^-l)(/^-7)]^/^ 



(44) 



The last equation thus gives the gap H23|l found in the 
random phase approximation. 
In the broken phase m ~ h and 



= - 



1-7 



+ [(i-/.2)(i-7; 



,1/2 



A(i^) = 2[(l-/.2)(l-7)] 



1/2 



(45) 
(46) 



In this case, one must be careful about the physical in- 
terpretation of A^^). It is not the true gap of the sys- 
tem, opened by the tunneling between the two classical 
ground states, which is knownffi to be exponentially small 
exp(— aA^)]. It is thus not the quantity denoted as A 
in Fig.^ but the one denoted as A', which is the excita- 
tion energy in the vicinity of one of the classical magne- 
tization [up to some terms of order exp(— aA) which are 
not accessible in the 1 /N expansion] . It is thus obtained 
by computing the energy difference between the ground 
state and the second excited state. For the same reasons, 
the ground-state energy as given by Scq is only valid up 
to exponentially small terms. 

Finally, we can compute the first quantum correction 
to the observables by using Eq. (|77|l and the expressions 
of the S operator truncated to lowest order and expressed 
with the &'s 



Sx ^ (flt + a) 



2 Vl-e 
iAi/2 



1/4 



(5t+6) 



^ - 2 

iAi/2 /l-e\i/4 



2 \l + e 



N I f 1 

= — + - 1 - 



(6t - b) 



2 2 V VT^/ 

e 



1 



-.b'^b 



2vr 



(47) 



(48) 



(49) 



The one spin averages in the rotated basis deduced from 
these formulas read 

2(^,)/A^ = 2(^y)/Ar = + O(l/A^), (50) 

2(^,)/Ar = l + l(^l--^=l=^+0(l/iv2). (51) 

In the symmetric phase, these expressions can be used 
by simply replacing the S by 5*, and by assigning the ex- 
pression H38|l to £. To obtain the results in the broken 
phase, one should pay attention to the fact that the cor- 
rections to {Sx) and {Sy) on the one hand and to (5*2) on 
the other, are not of the same order. One then gets 



2{Sx)/N 

2{Sy)/N 

'2{Sz)/N 



0, 

TO, 



(52) 
(53) 
(54) 



which are nothing but the mean-field results. In order 
to obtain the first quantum correction to these results 
in the broken phase, one must also be able to compute 
the 1/iV corrections to 2{Sx)/N and to 2(Sy)/N. Such 
a calculation requires that we go beyond the Bogoliubov 
transformation and will be presented in Sec. It is 
nonetheless possible to extract some of the first quan- 
tum corrections to spin-spin correlation function in both 
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phases 



A{Sl)/N^ 

aCsD/n^ 



1 ( i±£ 

TV VI 

1 fl-e 
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1/2 



1/2 
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r2 



4(5,5,) 7iV^ 
4(5,5,)77V2 ^ 

4(5,5,)77V2 = 0(1/7V). 



(55) 

(56) 
(57) 

(58) 

(59) 
(60) 



Here again, the BogoUubov transformation is a too simple 
calculation to obtain all the corrections. In the sym- 
metric phase however, we know that the spin-flip sym- 
metry is not broken, so that Eqs. H59|) and 16U|) in fact 
vanish. In the broken phase, although Eq. H6U|) still van- 
ishes, Eq. (|59|l is nonvanishing. We refer the reader to 
AppendixOwhere we have gathered the first 1/iV correc- 
tions (computed as explained in the next section) . Let us 
mention that sticking to the results (|55|l - (|6(J|I presently 
available, the only two-spin correlations we can compute 
in both phases, and better than in the mean-field treat- 
ment is 



4{Sl)/N' 



1 fl-e 



1/2 



(61) 



Luckily, as explained in Sec. IVII Al this is precisely 
the correlation function that is required to compute the 
rescaled concurrence in the parameter regime ft, > for 
7 < or /i > ^ for 7 > [see Eq. (I1U5|| ]. However, 
when 7 > and h > ^/7, one needs to know the 1/N 
corrections to (5,) and (5^) [see Eq. Hl()8(l ]. which can 
not be computed with the Bogoliubov transform. 

Of course, it would be tempting to go beyond these 
first-order corrections pushing the expansion of the spin 
operator to the next order. Unfortunately, from the or- 
der (l/iV)i/2, the Hamiltonia n is no more quadratic and, 
consequently, cannot be diagonalized by a Bogoliubov 
transformation. This clearly calls for a more sophisti- 
cated method which constitutes the aim of this article 
and is presented in the next section. 



V. HIGHER ORDER QUANTUM 
CORRECTIONS 

In this section, we show how to to go beyond the first- 
order quantum corrections given previously. We explain 
how to compute higher-order corrections in the ex- 
pansion thanks to the CUTs technique. These results are 
used, together with a scaling argument, in order to ex- 
tract the finite-size scaling exponents of various physical 
quantities at critical coupling. The analytical calculation 
of these exponents is the central result of this study. 



Introduction to continuous unitary 
transformations 



To compute higher order terms in the expansion, 
we need to diagonalize the Hamiltonian H28|) truncated 
at the corresponding order, i.e. to perform a unitary 
transformation that brings the Hamiltonian to a diag- 
onal form. Finding such an unitary transformation is 
clearly not an easy task. It was the idea of Wegner— ten 
years ago, and independently of Glazek and Wilson22i2S 
to replace the single step unitary transformation by an 
infinite product of infinitesimal transformations, that is 
by CUTs. In this way, one does not have to solve an 
algebraic problem but a set of differential equations. 

We shall now briefly recall the formalism of the CUTs, 
and refer the reader to Ref. \^ for a pedagogical intro- 
duction and more references on this topic. The idea of 
the CUTs is to diagonalize the Hamiltonian in a contin- 
uous way starting from the original (bare) Hamiltonian 
H = H{1 = 0). A flowing Hamiltonian is thus defined by 



H{1) = U\l)H{0)U{l), 



(62) 



where / is a scaling parameter such that H{1 — oo) is di- 
agonal and U{1) is a unitary transformation, i.e., satisfy- 
ing U(l)W{l) = W{l)U{l) = 1. A derivation of Eq. ^ 
with respect to / yields the differential equation (called 
the flow equation) 



diH{l)^[r^{l),H{l)l 



(63) 



where the generator of the unitary transformation ri{l) is 

T^[l) - diU\l)U{l) - -U^{l)diU{l). (64) 

Let us remark that the spirit is a bit different to what has 
been done with the Bogoliubov transform. Indeed here, 
H{oo) is unitary equivalent to (0) and is diagonal in the 
original basis in which H{Q) is nondiagonal. Conversely, 
when using the Bogoliubov transform, we expressed the 
same Hamiltonian H{Q) in two different basis (the a's 
and the &'s). 

One must also simultaneously keep track of the change 
of basis on other operators in which one is interested. De- 
noting by = ^{l ~ 0) such an operator in the original 
basis, one defines the flowing operator 



= [/^(/)f}(o)c/(0, 



(65) 



which is subject to the same flow equation as H , namely, 

dm)^[i^{i),n{i)]- (66) 

This allows one to compute the expectation value of 
any operator 11 on an eigenstate \'ip) of H. Indeed, one 
has 

{i^Mii) = {i}\U{l = c5o) = oo) U^l = oo)\i'), (67) 

where W{1 — oo)\ip) is simply the eigenstate of the diag- 
onal Hamiltonian H{1 ~ oo). 
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We still have to give a prescription for the anti- 
hermitian generator r]{l), that has to be chosen to bring 
the final Hamiltonian into a diagonal form. Wegner pro- 
posed 7?(0 = [ffd(0,ffod(0] = [H^{l),H{l)], where 
and Hod are the diagonal and off-diagonal part of the 
Hamiltonian. Other generators have been proposed, and 
we will use the so-called quasiparticle conserving gener- 
ator which is better adapted to our problem. We shall 
now give the motivation for this choice and then derive 
the flow equations. 



B. The flow equations 

When expressed in terms of the bosonic operators a 
and at, the initial Hamiltonian l|28|) reads 



H{0) = Ho{0) + H+{0) + H^{0) + H+{0) + 



H^{0), (68) 



where ^ — (^1^2)^ '^^^ '^■^ 1' '^^ ^ subscripts indicate 
the number of created (-f ) or annihilated (— ) excitations. 
When using Wegner's generator, one would take -ffo(0 as 
the diagonal part and the rest as the off-diagonal part so 
that the generator would be 




FIG. 5: Structure of the Hamiltonian in the bosonic operators 
a and a^, for terms involving less than six a or a) operators 
only (see the main text). 



^{^) = [Ho{l), {I) + {I) + (l) + H2 (01- (69) yet exphcit. Both points arc remedied writing 



Such a choice of generator suffers from one major draw- 
back: the flowing Hamiltonian H{1) docs not remain 
band diagonal as it is initially, but contains terms creat- 
ing or annihilating any number of excitations , k > 1. 

A possible choice of generator avoiding this problem 
is the so-called quasiparticle conserving generator first 
proposed by Mielke^ in the context of finite matrices 
and generalized to many-body physics by Knetter and 
Uhrigifi. For the problem at hand, this generator reads 



V{1) = H+{1)-H^{l)+H+{1) 



H-{1), 



(70) 



and coincides in the symmetric phase {Hf- = 0) with the 
generator proposed by Stein^. 

From a practical point of view, this generator is also 
more convenient than Wegner's generator because one 
does not have to compute a commutator to find the gen- 
erator. In addition, the flow equations are quadratic in- 
stead of cubic in the Hamiltonians 



dMl) = 2[[H+{l),H^{l)\ + [H+{l),H^{l)\)in) 
diH+{l) = [H+il),Ho{l)]+2[H+{l),H^{l)], (72) 
m+il) = [H+il),Hoil)]. (73) 

One can furthermore check, for example, that no three 
particle term is generated during the flow since the sum 
of [H+{l),H+{l)] and [H+ {I) , H+ (l)] judiciously cancels. 

Equations H7I|1 - H73() are unexploitable as they stand, 
because we can only follow the flow of coupling constants, 
not of operators, and because the 1/N expansion is not 



Hoil) 



]\fa+S-l ■ 



E 

2^ ]\fa+S-l/2 ' 
a,i5eIN 



E 



]\[a+S ■ 



(74) 
(75) 
(76) 



where Aa — a^^a". The couplings ^^'^(0 are precisely 
the flowing quantities we wish to follow. The fc = 0, 1, 2 
subscript is associated to the number of created excita- 

4. k 

tions a' . The a subscript keeps track of the number of 
and a operators in the operator. Finally, the i5 
superscript codes for the successive 1/iV corrections. For 
example, the diagonal part of the Hamiltonian is 

Fo(0 = iv<^(0 + 4!^(0 + ^C(0 + --- 



N 



C(0 



2 2 



a 



(77) 



where only terms up to order 1/N have been written 
explicitely. 

To make things even more concrete, let us write the 
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initial conditions {I = 0) 

C(0) = 2 " ^^^^ 

4^^(0) = '-^ = Se,^^\ (79) 

4'o(0) = 0, (80) 

4°1(0) = 2 + 2;im-3m2-7 = A(^)(oo), (81) 

4^1(0) = -2 + 3jn^+j, (82) 

4°2(0) = -2 + 3m2+7, (83) 

and all coefficients of higher order vanish. 

We have depicted the structure of the Hamiltonian in 
Fig. [SI for all terms involving less than six a or opera- 
tors. Let us emphasize that this figure does not represent 



The C coefficients arise from the computation of the com- 
mutators 

= EC:^>at^'"'"A„,+„.+,._„if/>/', (87) 

n 

E^i' .j''n A 3 "3 -c ■! ^ -n 

Ci',a" Aa'+a"+j'-„a if J <J , 

n 

and are equal to 

c£/{> = nl {Cl,C-„ - , (88) 

C2 being the usual binomial coefficient a\/[n\{a — n)!]. 
The various sums in (|84|) - (|86|l are constrained by the fact 
that all subscripts and superscripts have to be positive. 
Thus for example, in the first sum of Eq. H84|l . n runs 
from to 1 -I- (5, a' from to a — 1 -I- ?7, and 6' from to 
5+l-n. 

We shall now explain how to diagonalize the Hamilto- 
nian order by order in 1/N . The most important point is 
to ensure that corrections to a given order are unchanged 
when going to next orders. This means that terms of 
order (1/7V)~^ and (1/A^)~^/^ must not appear in the 
flow equations. Here, this is indeed the case because (i) 
the term of order (l/iV)~^ is proportional to the iden- 
tity operator and (ii) the term of order (1/A^)~^/^ has 



the Hamiltonian in the eigenbasis {|n)} of the number 
operator a^a. Let us explain on the example of the oper- 
ator a^- at a how it should be interpreted. This operator 
annihilates the states |0) and |1) containing zero or one 
excitation. The first state it does not annihilate is |2) 
and its action is to increase the number of excitations by 
one, and so gives a state proportional to |3). For these 
reasons, this operator is written on the first upper off- 
diagonal, at the position (3, 2). It is, however, clear that 
this many-body operator also acts nontrivially on any 
state \n>2). As another example, the zeroes shown in 
the figure mean that the Hamiltonian does not contain 
terms creating three or more excitations. 

Equations H74|l - H76() can be inserted in H71I73|I to give 
the flow equations of the various h!f^^{l) couplings 



(84) 
(85) 
(86) 



been suppressed by the appropriate choice of ^o- Conse- 
quently, couplings of order (1/A^)~^ and are 
simply /io°o(0 = eo(oo) and hfl{l) = for all /. 
At the next order (1/^)°, Eqs. (IHll)-® give 

dih^oiil) = (89) 
az<l(0 = -8/4f (0 (90) 
5,4°)(0 = -2/4°l(/)/4°)(0. (91) 

These equations are simply the Bogoliubov transforma- 
tion, written in a differential forn;^. It is a simple matter 
to recover Eqs. (|41|) and H42|) from the two constants of 

the flow, 2h^oM0 - 4°1(0 and 4a (0' " 4/i^°^(02. 

The same type of calculation can be performed for the 
observables Sx , Sy and Sz as well as the correlation func- 
tions SaSjS. As the analysis closely resembles the one we 
performed for the Hamiltonian, and as we want the main 
text to be as fluid as possible, we have gathered the flow 
equations of these observables and their derivations in 
Appendix ^ Let us simply mention that contrarily to 
the Hamiltonian, the observables do not retain the sim- 
ple structure they exhibit in the beginning of the flow, 
even when using the quasiparticle conserving generator. 

The flow equations for the Hamiltonian and the ob- 
servables can be integrated exactly in both phases, at 



I 

^ ^ '^a',a-a'-l+n"'l,a'\''!"'l.,a-a'-l+n\'') '^a' ,a-a' -2+n'^2,a' \'' )"'2,a-a' -2+71^ ) ' 

I 
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least up to the orders considered in this work. From the 
solutions, one can compute the renormalized quantities 
(i.e.; the values at / = oo). Note that the possibility to 
get exact corrections may originate from the inte- 
grability of the LMG model and is not a generic feature of 
the CUTs. We explain how to compute the solutions to 
the flows of the energies and obscrvablcs in Appendix IbI 
The results for the renormalized quantities arc gathered 
in Appendix lUl for the symmetric and broken phases. In 
the next subsection we explain how to make use of these 
results to compute finite-size scaling exponents. 



C. Computation of finite-size scaling exponents 

In this section, we explain how to make use of the 
results of the expansion to compute the finite-size 
exponents of the spectrum and of the correlation func- 
tions of the anisotropic model. Let us first focus on the 
results obtained in the symmetric phase, see Appendix 
EH iFromEqs. (ITTTOll . dTTTTll . ^2^, HCM . iHHHl) 

one sees that any physical quantities $ that we have com- 
puted can be written as 



$jv(/i,7) 



(92) 



where the superscripts "reg" and "sing" stands for regu- 
lar and singular, respectively. A nonsingular contribution 
is understood to be a function of h which is nonsingular 
at /i = 1, as well as all its derivatives. As an exam- 
ple, if $ = 2{S^)/N [see Eq. (fUTTj) ]. the regular part is 
1 -|- and the remaining forms the singular part. Now 
let us suppose h is close to its critical value 1. Then 
all terms involving the polynomial functions become 

small compared to the terms involving the P^'^^s, by a 
factor R(h,"f) = [h — l)(/i — 7). In this limit, one can 
write 



S(fe,7) 
TV"* 



7VS(/l,7)3/2,7 



(93) 



where J-^ is the scaling function for the physical quantity 
$, depending only on the scaling variable A^S(/i, 7)'^/^ 
[but not separately on N and S(/i,7)], ^$ and n$ are 
exponents that wc list in Tabic Now we use the fact 



$ 




n# 


n<f, + 2^$ /3 


eo 


1/2 


1 


4/3 


A 


1/2 





1/3 


2{S.>/iV 


-1/2 


1 


2/3 




-1/2 


1 


2/3 




1/2 


1 


4/3 




-1/2 


1 


2/3 



TABLE L Table of exponents (see the text for their significa- 
tion) for the physical quantities $ considered in this work. 



that there can be no singularity in any physical quan- 
tity in a finite-size system. This implies that the sin- 
gularity of S(/i,7)^* has to be canceled by the one of 
JF* [-/V5(/i, 7)'^/^, 7] . Thus one must have JF<i,[2:,7] ~ 
2;-2?*/3^ which in turn implies the following finite-size 
scaling 



$r'(/i = i,7) 



TY»<i.+2C4./3 ' 



(94) 



where the exponents n$ -|-2^$/3 are listed in the last col- 
umn of Tabled and the a$'s are constants which cannot 
be determined with this scaling argument. These finite- 
size exponents arc compatible with the values obtained 
from numerical diagonalization and discussed in the next 
section (see Figs. El and [Tj). 

Let us finally mention that the exponents for the 
ground-state energy per spin eg and for the gap A can 
also be computed following the same line of reason- 
ing, from the results obtained in the broken phase [see 
Eqs. HC39|) and (jC40|) ]. The only difference is that one 
has to replace the variable S(/i,7) — (h — l){h — 7) by 
^^(/i) = \ — h? . For the observables we did not go to suf- 
ficiently high order to perform the same analysis. How- 
ever, we can argue that from the results for eo and A that 
the scaling variable is A^4'(/i)'^/^. Then, one can easily 
show that we recover all exponents listed in Table |I| 



VI. NUMERICAL RESULTS 

In order to check the validity of the approach used here 
to compute the critical exponents, we have numerically 
diagonalizcd H in the sector 5* = N/2 where the low- 
energy states lie. The dimension of this subspace is + 1 
but the spin-flip symmetry (0J further divides by 2 the 
dimension of the matrices to be diagonalizcd. These sim- 
plifications allows one to easily investigate rather large 
system size (up to iV = 2^° spins for the energy spectrum 
and N = 2^^ spins for observables requiring eigenstates). 
For each system size N , we have computed the ground- 
state energy eg, the energy gap A, the magnetization 
{Sz) and the correlation functions (5*^), {Sy), and {S'D 
at the critical point h — \. For clarity, we only present 
the results for 7 = but we have checked that the scaling 
exponents are the same for any 7 ^ 1 as expected from 
our analysis, the peculiar case 7 = 1 being discussed in 
SecHTBl 

As can be seen in Fig. [S] the leading nontrivial finite- 
size corrections for eg and A are indeed proportional to 
N-^/^ and N'^/'^. Note that this latter result was al- 
ready observed numerically by Botct and JuUien^ twenty 
years ago but, to our knowledge, never proven by analyt- 
ical arguments. 

Our results are in contradiction with a recent study 
by Reslen and co-workersi^ who computed numerically 
some correlation functions. They indeed obtained that 
2{Sl)/N - 1 ~ iV-0-55±o oi. The source of this discrep- 
ancy comes from the too small system size they investi- 
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FIG. 6: Log2-log2 plot of the ground-state energy per spin eo 
and the enery gap A at the critical point ft = 1 as a function 
of N, for 7 = 0. 
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FIG. 7: Scaling of various observables as a function of A'' 
(logj — log2 plot) at the critical point h = 1, for 7 — 0.. 



gated (100 < N < 500). Indeed, in this range, we have 
checked that the scahng exponent is indeed close to 0.55, 
but for larger size as those considered here {N ~ S^''), 
we found numerically that the fitted exponent gets closer 
and closer from 2/3 which is the value predicted by our 
approach. 



VII. ENTANGLEMENT PROPERTIES 

The entanglement properties in the LMG model have 
recently attracted much attentioni^iSiiSii^SiS. For the 

ferromagnetic case considered here, several quantities 
have been used to characterize the ground-state entan- 
glement. First, the concurrence^ which measures the 
two-spin entanglement has been shown to develop a cusp 
at the critical point^ with a nontrivial scaling behav- 
ior that we shall discuss below. Secondly, the von Neu- 
mann entropy has also been computed as a function of 
both the anisotropy and the magnetic fieldSi. Its scaling 
at the critical point has suggested the possibility of an 
underlying conformal theory describing the LMG model 
with central charges c = 3/2for7 = l and c = 1 other- 
wise. Finally, we wish to clarify the situation about an- 
other measure proposed by Somma et alm^, the so-called 
SU(2)-purity. For a system where all spins are equivalent 
as in the LMG model, this quantity is nothing but the 
square of the reduced Bloch sphere radius defined by 



N 



(95) 



The latter relation is due to the spin-flip symmetry Q 
and is only valid for the exact finite- ground state. It is 
thus simply connected to the magnetization in the field 
direction so that the SU(2)-purity does not bring any new 
information about the entanglement of the ground state in 
this model. Furthermore, this purity is trivially computed 
in the thermodynamical limit using the semiclassical de- 
scription discussed in Sec. IIIII Unfortunately, for the 



microscopic degrees of freedom which are the spins half, 
there is no entanglement in this limit since the ground 
state is a pure product state [see ansatz (|19|l ]. Clearly, 
Somma et al. did not understand that the richness of 
the entanglement in this model was given by the finite- 
corrections which we shall now discuss in the framework 
of the concurrence. 

Introduced by Wootterai^, the concurrence which mea- 
sures the entanglement between two spins half is defined 
as follows. For a given pure state we define the den- 
sity matrix p = \ip){il)\ and the two-spin reduced density 



matrix obtained by tracing out p over all spins ex- 
cept spins i and j. Next, we introduce the spin-flipped 
density matrix pi^j = {ay ® ay) p*j (ay (g) ay) where p*^ 
is the complex conjugate of pi_j. The concurrence is then 
defined by 



Cij = max {0, pi~ P2~ ^4} : 



(96) 



where the pk 's are the square roots of the four real eigen- 
values of Pi J Pi J, and where pi > pi+i- The concurrence 
vanishes if the reduced density matrix p^.j can be decom- 
posed into uncntangled pure states whereas it reaches 
its maximal values 1 for maximally entangled (two-spin) 
states such as the famous Einstein-Podolsky-Rosen state 



|EPR) = -i=(|Ti) 



lit) 



(97) 



In the LMG model, since all spins are equivalent, the 
concurrence does not explicitcly depend on i and j and 
we shall omit these indexes in the following. In addition, 
since the ground state lies in the sector S = N/2 the 
concurrence can be easily expressed in terms of the ob- 
servables computed in the previous section. Indeed, for 
any symmetric state, the reduced density matrix reads^: 



u 



w+ .tZ|_ a;;;}. 

x+ y y X* 

x+ y y X* 

u x_ x_ u_ 



(98) 
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with 



v± 
x± 

y 

u 



N^-2N + 4(5,^) ± 4{N - 1){S,) 

4:N{N- 1) 
{N~1){S+)±{[S+,S.]+) 
2N{N-1) 

AN{N 
N{N -1)' 



(99) 
(100) 
(101) 
(102) 



For an eigenstate of the spin- flip operator as, for instance, 
the ground state of the LMG model at finite N, one fur- 
ther has x±=0, so that the concurrence is simply given 



C 



2 max(0, \u\ — y) if 2y < ^v+v- + \u\, 

2 max(0, y — ^v+v-) if 2y > ^/vyvZ + \u\. 



since for I7I < 1, one has 1(5^ — Sy)\ > 0. It is worth 
noting that, obviously, the concurrence C goes to zero 
in the thermodynamic limit since the ground state be- 
comes a product state. However, the relevant (nontrivial) 
quantity for this problem is the rescaled concurrence Cr,. 
This rescaling takes into account the fact that the en- 
tanglement is equally shared between all spins and shall 
be interpreted as the connectivity of each spin. Using 
expression H61|l obtained by the Bogoliubov transform, 
one can then compute the asymptotic behavior of the 
concurrence in both phases. One thus has 



(103) 




(106) 
(107) 



A. The thermodynamical limit 

Using the expressions of the observablcs gathered in 
Sees. IC II and IC 21 one can show that for ^ < h, one 

has 2y < ^v+v- + \u\ so that the rescaled concurrence 
Cr = (iV - 1)C reads^ 

^H^^' = |( K^' - I - + (5f ) ) , (104) 
which can be rewritten as 

Cf^' = 1 - (105) 



These expressions are valid for any I7I < 1 and coincides, 
for 7 = 0, with the results found by Reslen et ali^. They 
also prove that a real cusplike singularity occurs at the 
critical point as inferred from numerical results in Ref^. 
Note that in the broken phase, the concurrence vanishes 
for h = ^/7 which is precisely the point where the ground 
state is a pure product state at finite N (see Sec. Illljl . 
At this special point one has 2y ~ ^/vyuZ + \u\ and for 
h < .^7, one finds that 2y > ^v+v- + \u\ so the second 
expression in Eq. H103|) must be considered to compute 
the concurrence, namely 



N{N -2) + 4{S^,) - 4{N~1){S,) 



(108) 



Using the expansion of {Sz 
one then finds: 



icTTll and (S'2) dcTTIl . 



1 



1-7 

l-/l2- 



(109) 



We emphasize that for ft, = 0, this latter expression shows 
that the rescaled concurrence of the ground state is al- 
ready nontrivial since it only vanishes for 7 = 0. Further- 
more, let us stress that this result is beyond the reach of 
the Bogoliubov transform. The results for the rescaled 
concurrence in the thermodynamical limit are summa- 
rized in Fig. |Hlfor an anisotropy parameter 7 = 1/2. 

It is important to note that, in the broken phase, we 
have used the expressions of the correlation functions 



r 



computed in a 1/N expansion around the ansatz state 
\'4){9[),(f)Q)) [see Eq. (I19|l ] which is not the exact ground 
state since it has, in particular, nonvanishing values of 
{Sx) and {Sy). The exact ground state jV'o) and the ex- 
act first excited state j^/'i) indeed belong to subspaces of 
fixed and different spin-flip symmetry. However, in the 
broken phase, correlation functions for these states are 
identical up to exponentially small terms ~ exp(— oTV). 
Thus, for any state Y^a) = cosajT/io) + sinal-i/"!), the av- 
erage of an operator Vl such that [il, cr*] = (as, for 
instance, S^, Sz or S'^) is independent of a and satisfies 



{■^{a)\n\^{a)) = (^olf^lV-o) = (^ilf^lV-i). (110) 
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FIG. 8: Rescaled concurrence as a function of the magnetic 
field, in the thermodynamical limit, and for 7 = 1/2. 

This identity is valid up to exponentially small terms 
and justifies our use of correlation functions computed 
from the 1/A^ expansion around one of the two broken 
symmetry state. 




0.5 1 1.5 



h 

FIG. 9: Rescaled concurrence as a function of the magnetic 
field, for 7 = 0, and for system sizes iV =16, 32, 64, 128, 256 
and 00 (bottom to top). Arrows indicate the behavior of the 
finite size correction in various regions. 

Goldstone mode associate to the broken rotational sym- 
metry along the z axis. 



B. Finite-A^ corrections 

Let us now discuss the finite-A^ corrections to the 
rescaled concurrence. For I7I < 1, apart from the critical 
point, the leading corrections are proportionnal to 1/A^ 
since the expansion of observables is not singular. A con- 
trario, at h = 1, using Eq. H105|) and results of Sec. lC 1 el 
one gets 

C't' - I - ayyN-'/^ (111) 

which is indeed the result suggested in Refsii^i^i. We 
have summarized these scalings in Fig. Elfor 7 = 0. 

C. The isotropic case 

For completeness, let us briefly discuss the isotropic 
case. As explained in Sec. IIIBI for 7 = 1, the ground 
state is simply the Dicke state {N/2, Mo)} whose rescaled 
concurrence is 

^{N^-UI§)[{N-2y-4Mi]y 

In the thermodynamical limit, one thus has a discontinu- 
ity at the critical point where the rescaled concurrence 
jumps from (for /i > 1) to 2 for {h < 1), and goes to 
1 in the zero field limit. Let us finally mention that this 
behavior can not be recovered by taking the limit 7 — > 1 
of the anisotropic results. Indeed, although Eq. 1] 106(1 
properly predicts Cr = in the symmetric phase, the 
broken phase result H1U7|) diverges in the limit 7 — > 1. 
The physical origin of this divergence is the presence of a 



VIII. CONCLUSION 

The aim of this work was mainly twofold. We have 
computed the finite-size scaling exponents for the spec- 
trum as well as one- and two-spin correlation functions 
of the quantum critical LMG model. We have also 
characterized the two-spin entanglement properties of its 
ground state by computing the concurrence. 

Such studies had already been performed numerically, 
by computing the exact eigenenergies and eigenvectors. 
The numerical calculations are not too demanding since 
the size of the matrix to be diagonalized grows linearly 
with the number of spins N. Such numerical results were 
presented in Sec. IVII where we computed the eigenener- 
gies for system sizes as large as A^ = 2^°, and the corre- 
lation functions for less than A^ = 2^^ spins. 

Our main goal however was to provide an analytical 
derivation of the simple finite-size scaling exponents ob- 
served numerically (Sec. IV Cp . as well as formulas for 
the rescaled concurrence in the thermodynamical limit 
(Sec. IVII A(l . Both of these calculations require the 
knowledge of 1/A^ corrections to the thermodynamical 
limit (mean-field) solution (Sec. Illljl . On the one hand, 
as far as entanglement properties are concerned, it is clear 
that one has to go beyond the mean-field wave function, 
which is separable and can thus display no entanglement. 
The nontrivial entanglement properties arise in the finite- 
size system, and the physically interesting property is the 
first 1/A^ correction to the concurrence. On the other 
hand, finite-size scaling exponents at criticality can be 
extracted from the knowledge of the singularities aris- 
ing in the 1/N expansion, by making use of a finite-size 
scaling argument. 

In Sec. CH it was shown that the Bogoliubov trans- 
form solution of the bosonic Holstein-Primakoff represen- 
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tation of the Hamiltonian truncated to quadratic order 
does not capture enough of the 1/A^ corrections needed 
in this program. One thus has to work with a non- 
quadratic Hamiltonian which, in turns, means that one 
needs a method that can solve such interacting Hamilto- 
nians. We have used what we think is the most efficient 
and elegant method to perform this task, namely the 
CUTs fSec. IV A)| . With this method, the spectrum and 
correlation functions are the infinite time limit of some 
renormalization flows. For the LMG model, such flows 
can be computed exactly, yielding exact successive 1/N 
corrections. The integrability of the flows, which is not 
a generic feature of the CUTs, may originate in the fact 
that the model is solvable by algebraic Bethe ansatz. One 
should however stress that it is a difficult task to com- 
pute 1/iV corrections for the spectrum from the Bethe 
ansatz solution2^i2&, and even more complicated to do it 
for correlation functions. 

The combination of the 1 /N expansion in the Holstein- 
Primakoff representation, the CUTs and a scaling argu- 
ment, opens the possibility to compute finite-size scaling 
exponents as well as some nontrivial 1 /N effects for other 
models as well. Of particular interest is the reduced BCS 
model used to describe superconductivity in ultrasmall 
grainsi^ This will be the object of a forthcoming publi- 
cation. As they stand our results are also directly appli- 
cable to the Dicke model, which in the zero-temperature 
limit can be put in a one-to-one correspondence with the 
LMG modelii^ We hope more models will turn out to be 
analyzable by following the same lines of reasoning. 
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APPENDIX A: FLOW EQUATIONS OF THE 
OBSERVABLES 



Let us denote by a general flowing hermitian ob- 
servable. This observable can be decomposed in the fol- 
lowing form: 



n{i) = n^{i) + Y,Mii) + ^'k{i)\ 



(Al) 



where the sum runs over all nonnegative integers /c's and 

fi^T = ^- The flow equations are found by computing 
[ri[l),n{l)]. One finds 



din2{i) 



[H+{l),n^{l)]+h.c.j (A2) 

.[[H+{i),n^{i)] +h.c.), 

ff+(0,f^o(0] + [H+{l),n^{l)]^ (A3) 

H+{i),nt{i)] + [H+{i),n-{i)]^ (A4) 

H+ii),nUii)] + [if+(0,O-+i(/)]^A5) 



In a notation similar to Eqs. H74() - (|76() . we have intro- 
duced 



Ml) 



E 
E 



J\fa+S-l ' 



(A6) 
(A7) 



Inserting these formulas into Eqs. l|A2p - ljA5p yields 



(i5-i5' + l-n) 



("5') ^,v,,(<5-<5' + l-») ' 

2,Q — a' — 



(0 



71. a.' ,8' rL,a' ,8' 

,("5') I'M, + , ^2,3,n AS' ) , (5-S' + l-n) * 



(A9) 



/ . ^a',Q-a'-l-|-n"2,Q'V''7^1,Q-a'-l+nW "T '^a\a-a' -2+11 2,a'\'' >^Z,a-a' -2+n V'^ 

n^a',S' 71. a' ^5' 

E '^:^:'+n(r)h^2'^o^-':^^tw+ E (aio) 



E^2,fc-2,n (,("5') ^n, ,(<5-<5' + l-'i) I ^2,fc+2,ri AS'),,. (S-S'+l-n) ^ , 

'^a',a-a'+n"-2,a'\^)^k-2,a-a'+n\'-) ^ a' ,a-a' -2+n"-2,a'\'' )^k+2,a~a' -2+n \^ ) 



(<5') (<5-<5' + l-n) 



n^a'S' 



In the above equations, the asterix denotes complex con- 



jugation, the real part and we have defined the coeffi- 
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cients 

= «■! (C^"' - C2„ ) , (Al 1) 

arising from the calculation of the commutator 

n 

(A12) 



The initial conditions for Sz^ for example, are 



- 0) = 1/2, <i(/ = o) 



(0)/ 



-1 



(A13) 



and all other coefficients are zero. Once the solution of 
the flow equations for the Hamiltonian are known, the 
flow equations for the observables can be solved. 



APPENDIX B: INTEGRATION OF THE FLOW 
EQUATIONS 

In this appendix, we explain how to solve the flow 
equations, order by order in the 1/A^ expansion. For 
the observables, we will focus on 5*2, since the computa- 
tions for Sx and Sy differ only with respect to the initial 
conditions (and also the solutions for Sz are found to 
be simpler, which makes this Appendix more readable). 
Furthermore, Sz{l = 0) has real coefficients, so this will 
remain true for any Z, since also has real coefficients. 
This allows us to drop the real part 3? in Eqs. HA8p - HA10|l . 
The same applies to S^, but not to Sy whose coefficients 
are purely imaginary. 



1. Order (l/iV)"\ {l/N)-'^/'\ and (1/7V)'' 

a. Flow of the Hamiltonian 

The two first orders for the Hamiltonian have already 
been discussed in the main text: h'^Q{l) = eo(oo) to order 
(l/N)-^ and h^°l{l) = to order {l/N^^^'^ Further, we 

recall how crucial it is to ensure that hflil = 0) = 0, so 
that this remains true at any /. In the following, we will 
thus drop all contributions to the flow equations involving 
this term. 

To order (l/A^)", the flow equations are given by 
Eqs. (IH3)-linil)- The equations for h'^olU) and h^^oi^) are 



easily solved by noticing that h'^olil)^ — 4:h'^o{l)^ is a con- 
stant of the flow. One gets the hyperbolic solutions 



4°o(0 = 



tanh [2AW(Z + io)] 

-sgn(£)A(R) 
2sinh[2A(R)(Z + /o)] 



A(B)/i(0, (Bl) 

= r(B)/2(0, (B2) 



where Iq is such that the initial conditions are fulfilled, 
namely A^^' = A^^V tanh [2A(^);o] . The quantities 



A<-^\ r(B), A<-^\ and e are defined in Sec.|TV| The func- 
tions ,fi{l) and /2(0 are defined by the two equations 
above, and satisfy /i(0) = /2(0) = 1. Finally, h'^l{l) is 
found thanks to the second constant of the flow, namely 

2<o(0-<l(0- 

There is also another useful representation of the so- 
lutions (jBip and ljB2p . Let us introduce a new "time" 
scale 



t = sgn(£)cxp \2A^^Hl + lo)] , 



(B3) 



with initial conditions now given at 

to =sgn(£)exp(2AW;o) • (B4) 
After some algebra to can also be shown to be equal to 

to = - f 1 + Vl-£^) ■ (B5) 



(B6) 
(B7) 



Equations IjBlj) and ljB2p now read 



t" + l 

t2~T' 



<l(t) A(«' 



The renormalized values at ^ — *■ oo are now found by 
taking the limit t — > too = sgn(e)oo. Let us remark that 
the off-diagonal coupling h^2oi^) &oes to zero and behaves 



as t for t 



This will be true for all off-diagonal 



couplings creating two excitations since the energy cost 
of such excitations, in the thermodynamic limit and for 
large t, is nothing but 2A^^\ so that the couplings must 
vanish as exp 



b. Flow of the Sz observable 

To order (1/A^)"^ the flow is 

di4°^{l) = Q, (B8) 

As a consequence, t^o°o(0 = '^o°o(0) = 1/2- 

To order (1/iV)^^/^, there is only one flow equation, 
namely 



(B9) 



Since for Sz the initial condition is wj°g(0) = 0, one gets 

Working to order {1/N)'^ one has three flow equations 
to solve, which are 



d^4]lil) ^ -2;i(>("^(0-44°)(0c.^°^(0, (BIO) 
diJ^lil) = -4;i(>f^(0-84°)(0a;^°)(0, (BU) 
d^J^lil) = /.(>("^(0-24°)(04°](0- (B12) 
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For the Sz observable considered here for illustration, 
the previous order solution tells us that the first terms in 
the three right-hand sides (RHSs) of the above equations 
vanish, so that the equations are homogeneous. Let us 
stress that it is not the case for Sx ■ One first has to solve 
the subsystem (|B11|I and l|B12(l . and then the RHS of 
(|B10|) is known so that it can also be integrated. In the 
t parametrization, one has 



^2^0 (too) 



ti-i 



ti-l 



(B20) 
(B21) 



Up to a factor 1/2, Eq. (|BT9|) is nothing but the 1/A^ 
term in Eq. (|5T1) . 



which is integrated in 



1 -t2 



1 -t2 



1-^2 



2+2 



4>) 



jt - to)il - tpt) 



(B13) 
(B14) 
(B15) 

(B16) 
(B17) 
(B18) 



satisfying the initial conditions a;Q^Q(to) ~ 0, ujl^l(to) = 

-1 and Lj!^°l{to) = 0. Fr om these, one can deduce the 
renormalizcd values 



1 



^2 - 1 



1 



, (B19) 



2. Flow of the Hamiltonian at order (1/A'^) 



The flow equations read 



1/2 



Cl(o+4/i^"^(o]/4!o(o 

hZil)+8h^2ii^)]hf}{l), (B23) 



with/i^°|(0) = 2h[^l{0) = 2m Vl - m^. The second equa- 
tion is independent of the first one so it can be solved 
separately. Then its solution can be inserted in the first 
equation. This first equation is then solved in two steps, 
the first being to solve the homogeneous equation, and 
the second being to find a particular solution of the in- 
homogeneous solution. Again, the simplest way to solve 
these equations is to use the t timescale. We simply give 
the solutions, which read 



i\/r~~ 



;t^/^t - l)l/2(to + 1)3/2(1 + _ 3i _ i^i) 



2m VT" 



m^ 



rtl/2(t-l)3/2fa + 1)5/2 
4/2(t0-l)3/2(t+ 1)5/2 



(B24) 
(B25) 



We see that both these off-diagonal couplings go to zero 
and behave as t~^/2 fQj- 1 too- This is a general feature 
of all off-diagonal couplings creating one excitation since 
their energy cost, in the thermodynamic limit and for 
large t, is nothing but A^^\ meaning that these couplings 
have to go to zero as exp(— A*^^^/). 



3. Next orders 

/^From order {1/Ny for the Hamiltonian and {l/NY^'^ 
for the observables, things really become more involved, 
and it does not make any sense any more to give the 
solutions in this Appendix. However, we will explain 
how one can find these solutions. 



Stein gave some explicit solutions for the flow of the 
Hamiltonian in the symmetric phase and in the I time- 
scale language [up to order (1/7V)^)^. The solutions 
he gives arc found to be polynomials in three functions: 
fi{l), /2(0 [see Eqs. (|BlJ and ^] and f3{l) = l + l [sat- 
isfying the same initial condition as /i and /2, namely, 
/3(0) = 1]. The degree of the polynom in /i can be re- 
stricted to be one, since one has (/) — e^/KO = 1 ~ 
In Stein's solutions, the exponential decrease of all off- 
diagonal couplings creating two excitations is found to 
mean that one can factor out a term /2(0- 

We have checked that to the maximal orders at which 
we worked, an ansatz of the following form: 

= ^o^'2[/i(0,/2(0,/3(0], (B26) 
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C(0 = /2^^(0n«[/i(0,/2(0,/3(0], (B27) 

C(0 = /2(0^S[/l(0,/2(/),/3(0], (B28) 

is always a solution provided the degrees in /2 and /a 
of the polynoms Pqi -Pii and P2 are large enough. The 
prefactor /2^^(0 in the solution for hi ensures the expo- 
nential decrease of these terms creating one excitation, 
with a time-scale being half the one for terms creating 
two excitations. 

In the t time-scale language, f2{l) reads 



rl/2. 



The P polynomial functions translate in fractions having 
the following properties. The denominator is of the form 
{t — 1)"" {t — 1)"+, with n- and n+ being integers. The 
numerator is a polynomial function in t and in ln(t/to) 
(this logarithm is /a up to multiplicative constant). 



/2(0 



(B29) 
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APPENDIX C: COEFFICIENTS OF THE 1/iV EXPANSION 



In this appendix, we give the first terms of the 1/A'^ expansion for various quantities. Note that we could have only 
given {H)/N = eo(A^), {Sz), and {SD which are sufficient to compute all these quantities (except the gap). Indeed, 
one has 

{H) = -'^{{SD + l{Sl)) 2h{Sz) + ^(1 + 7) (CI) 

and 

(s^) = {sD + {sD + {sD = ^ (^^ + 1^ . (C2) 

However, we have considered that explicit forms of (S*^), {Sy) are of interest because of subtle cancellations of some 
coefficients. 

For the special case 7 = — 1, we emphasize that CUTs have already been used by several groups^^iSiSi, notably to 
compute the 1/N corrections^. In addition, these corrections have also been obtained by a perturbative method^. 



1. Symmetric phase 



In what follows, we denote by ^(/i, 7) = (/i — l)(/i — 7). 



a. Ground-state energy per spin eo{N) 



eo{N) 



1 

iV2 



Fi^)(fe,7) Qi'\h,^) 
Eih,j) "^S(/i, 7)1/2 



+ (1-7)^1^ 



7V3 



1 

iV4 



_S(/i, 7)5/2 s(/i,7)2 
where the P's and Q's are polynomial functions of h and 7 that we list below 

P^'Hh,l) = i(l+7)(/i'+7)-^/i(l + 147 + 7'), 



24:h^ + (1 - 427 + + I67 



l(3) 



S(/i,7)4 ^S(/i, 7)7/2 



O 



N5 



(C3) 



(C4) 
(C5) 
(C6) 
(C7) 
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P^^\h,-f) = — (1 + 7) [32/1^ + (1 + 147 + 72)/i'' + 7(1 -667 + 7^)/i^ + 167= 



128 



20(1 + 3O7 + "f^)h^ + (1 + 367 - 8427^ + 867^ + j^)h'^ - 47^(11 - 547 + 117^) 



512 

-^(1 + j)h[96h^ + (1 - 347 + j')h^ 
1 

64 



647" 



(4 - 1847 + Aj^)h'^ - 7(5 - 2027 + 57^)/i^ - 167^ 



(C8) 



(C9) 



b. Gap A{N) 



with 



AiN) = 2S(/i,7)^/^ + ^ 



+ (1-7)' 



iV2 



S(/i,7) ^S(/i, 7)1/2 

Pf (fe,7) , Q^l\h,^) 
S(/i, 7)5/2^ S(/j,7)2 



1 



S(/l,7)4 ^S(/i, 7)7/2 



f^(^7) 



(^7) 



Qi'^(^7) 



^f^(^7) 



4'^(/^,7) 



2(l + 7)(/i'+7)-^/^(l + 147 + 7'), 



Q«(/i,7) = -(1+7)/^ + 27, 

i(l + 7)/i(/i^ - 7) - 18/i'' + (1 - 367 + j^)h^ + I67' 



— (1 + 7) \lUh^ + (7 - 627 + 772)/i'' + 7(7 - 2227 + 7-f^)h^ + 1127^ 
64 L 

-^(1 + j)h\^84h^ + (1 - 347 + -f^)h^ - 527^ 



O 



7V4 



25607/1* + (7 + 2127 - 45347^ + 2127^ + 7j^)h^ - 1287^(1 - I47 + 7^) 



1 r 

16 



(5 - 3147 + 5-/^)h'^ - 77(1 - 5O7 + j^)h'^ - 327^ 



(CIO) 



(Cll) 
(C12) 
(C13) 

(C14) 

(C15) 

(C16) 



c. One-spm expectation value (5z 



2{S, 
N 



with 



1 

N 



>P(fc,7) 
S(/i, 7)1/2 



S(/j,7)2 +S(/i, 7)3/2 



(2), 



1 



1 . 



1 

128' 



1 

7V3 



S(/i,7)V2+ s(/i,7)3 



(1 + 7) 56/i* - (1 - 987 + 72)/i2 



(3), 



(C17) 



(C18) 
(C19) 
(C20) 
(C21) 
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2Ah^ + (1 - 2267 + 7^)/i^ - 87(1 - 197 + 7^) 
3 



(C22) 



d. Two-spin expectation value (S^) 



with 



)(3) 

' XX 



= (/^-7) 
1 

+ 7V3 



1 1 



TV S(/i, 7)1/2 ^2 

pix^(fe,7) gi'l^(fe,7) 

S(/i,7)V2 + S(/.,7)3 



S(ft,7)2 +S(/i, 7)3/2 



4/1^ - 2(1 + ^"f)h^ + (1 + 87 + 372)/! - 47 



1 

2 
1 

64 



2/1^- (l + 37)/i + 27 



96(1 - 7)/i^ - 1287(1 - 7)/i^ + 4(3 - 657 + 897^ - 27-f^)h^ 



(C23) 



-(3 - 767 - 627^ + 1487^ - l?,"f^)h^ - 47(3 - 337 + 497^ - l9-f'^)h - 167^(3 - 27 - 7^) 
4(1 - l)h\2h^ - 27/1^ - 7(3 - 7)/i + 7(1 + 7) 



(C24) 
(C25) 
(C26) 

(C27) 



e. Two-spin expectation value {Sy) 



with 



iV2 



1 



S(/i, 7)5/2+ s(;,^^)2 



-7) 

O 



Q|/'J(fe,7) 
S(/i, 7)1/2 

1 

iV4 



^'i'^(/i,7) = 
Pi'n/i,7) = 



Ah^ - 2(5 + 7)/i2 + (3 + 87 + 7^)/! - 47^ 



2/1^ - (3 + 7)/i + 27 

1 r 



96(1 - 7)/i^ - 128(1 - 7)/i^ + 4(27 - 897 + 657^ - '6-f^)h^ 



(C28) 



64 

-(13 - 1487 + 627^ + 767^ - 37"^)/!^ - 47(19 - 497 + 337^ - 37^)/i - 167^(1 + 27 - 37^) 



4(l-7)/i 



2h^ ^2h^ + {l--i-f)h + -f{l + ^) 



(C29) 
(C30) 
(C31) 

(C32) 



/. Two-spin expectation value {Sl) 



7V2 



1 



Pii^(^7) 
S(/i, 7)1/2 



+ 2 



1 

iV2 



PL'Hfe,7) Qil^(fe,7) 

S(/l,7)2 ^S(/i, 7)3/2 



(2)/ 



(C33) 
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with 



P«(/i,7) = (l+7)-2/j = 2Pii)(/.,7), 



-^(l+7)/i 16/i2 + (l + 37)(3 + 7) 



(1 + 7)(3/i^ + l)-^h ih'^ + 3(1 + 7)2 



1 r 
2 



4/1^^ + (7 + IO7 + 7-f^)h'^ + AY 



64 



(1 + 7) 224/i^ + (13 - 2027 + rsY)h^ + 16Y 



-h 



32 
3 



112/1"^ + (59 - 1027 + 59-f^)h^ - 27(19 - 67 + 197^) 
3 



(C34) 
(C35) 

(C36) 

(C37) 

(C38) 



2. Broken phase 



In what foUows, wc wiU denote '^{h) = 1 — h^. 



a. Ground-state energy per spin eo{N) 



eo(iV) 



1 

iV2 



7/1^-2 + 7 (1-7)1/2 



2*(/i) 



1'(/l)l/2 



(C39) 



6. Gap A{N) 



A(iV) = 2(1 - 7)i/2^(/i)i/2 + i. 



-2 



(1 - 27)/i2 + 2 - 7 2(1-7)1/2 



1'(/l)l/2 



(C40) 



c. One-spin expectation values 



2(5.) 

iV 
2(5,) 

iV 

2(5,) 
iV 



*(/l)l/2 



AT 



(1 - 27)/i2 + 2-7 
2(l-7)i/2^(/i) + 5-(/i)i/2 



+ 



iV2 



(1-7)1/2/1 



*(/i)l/2 



*(/i)l/2 



iV2 



(C41) 
(C42) 
(C43) 



d. Two-spin expectation values 



4(5g; 

iV2 



fe2(;^2_3^27) 

(l-7)l/2^r(/i)l/2 



(/l2 -2 + 7)^(fe)l/2 

(1-7)1/2 



(C44) 
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4(5|) 
7V2 



1/2 



1 ^'(/l) 

iV(l -7)1/2 



1 

7V2 



2(1-7)1/^/1^ 
«'(/i)i/2 



+ (l_^)l/2^(/j)l/2 



(C45) 
(C46) 
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